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ORIGINAL ARTICLE 

Pharmacometric Approaches to Guide Dose Selection 
of the Novel GPR40 Agonist TAK-875 in Subjects With 
Type 2 Diabetes Mellitus 

H Naik 1 , J Lu 2 , C Cao 1 , M Pfister 2 , M Vakilynejad 1 and E Leifke 1 

The G-protein-coupled receptor 40 agonist (GPR40) TAK-875 is being developed as an adjunct to diet and exercise to improve 
glycemic control in patients with type 2 diabetes mellitus. Pharmacometric approaches such as model-based exposure- 
response and meta-analyses were applied to (i) characterize exposure/dose-efficacy responses of TAK-875, (ii) characterize 
the time course of glycosylated hemoglobin A1c (HbA1c) response with TAK-875 6.25 to 200 mg q.d. doses for 12 weeks, 

(iii) project and compare HbA1c response with dipeptidyl peptidase 4 (DPP-4) inhibitors and TAK-875 up to 24 weeks, and 

(iv) provide a quantitative rationale for dose selection in phase 3. On the basis of phase 2 data, relationships between TAK-875 
concentrations and HbA1c were well characterized by exposure-response models. EC 50 and E max of TAK-875 were estimated to 
be 3.16 |jg/ml and 0.366, respectively. Model-based simulations over 24 weeks indicated that the 25- and 50-mg q.d. doses of 
TAK-875 achieve efficacy as comparable with or better than that of commonly used antidiabetic agents. 
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Type 2 diabetes mellitus (T2DM) accounts for 90-95% of all 
diabetes cases; its global incidence is expected to increase 
from 171 million to 366 million by the year 2030. 12 Despite 
new classes of antidiabetic medications, the National Health 
and Nutrition Examination Survey has concluded that the 
proportion of patients who achieve a glycosylated hemo- 
globin A1c (HbA1c) level of <7% (the American Diabetes 
Association goal) remains relatively unchanged at about 
37%. 3 < 4 

TAK-875 (Supplementary Figure S1 online) though being 
developed as an adjunct to diet and exercise to improve gly- 
cemic control in patients with T2DM, differs from sulfonylu- 
rea-receptor activators and other secretagogues in its effects 
on p cells. Sulfonylureas and glinides stimulate insulin secre- 
tion even at low blood glucose concentrations; nonclinical 
data suggest that TAK-875 stimulates insulin secretion only 
at elevated blood glucose levels. 5 TAK-875 is an agonist of 
free fatty acid receptor 1 (or G-protein-coupled receptor 40 
agonist, GPR40) and is expressed in pancreatic islets cells. 
Binding of free fatty acids to the GPR40 receptor leads to 
glucose-stimulated insulin secretion. 56 This novel mecha- 
nism of action may allow safe and effective T2DM treatment 
by selectively improving glucose-dependent insulin secretion 
with a low probability of hypoglycemia. 

TAK-875 pharmacokinetics (PK) appear similar between 
healthy individuals 78 and patients with T2DM. 9 After single 
and multiple doses (25-800 mg) in healthy subjects and 
patients with T2DM, TAK-875 was absorbed from the gas- 
trointestinal tract with a median time to maximum con- 
centration of -3-4 h, and was slowly eliminated from the 
systemic circulation, with a terminal half-life of -28-51 



h. A dose-proportional increase in exposure (maximum 
concentration and area under the curve) without induc- 
tion of hypoglycemia was observed in placebo-controlled 
ascending single- or multiple-dose studies. 7-9 TAK-875 is 
primarily eliminated by hepatic metabolism via formation of 
glucuronide conjugates, with minimal renal clearance (CLr 
< 0.003 l/h). 7 

This study presents two pharmacometric approaches. First, 
utilizing phase 2 data, population exposure-response models 
were developed to (i) characterize the time course of TAK-875 
plasma concentrations and fasting plasma glucose (FPG) 
and HbA1c levels in patients with T2DM receiving 6.25-200 
mg oral q.d. doses of TAK-875, (ii) characterize the exposure- 
efficacy response of TAK-875, (iii) identify the sources of 
variability associated with PK and efficacy parameters of TAK- 
875, and (iv) project HbA1c response with TAK-875 up to 24 
weeks. Second, using publicly available clinical efficacy data, 
model-based meta-analyses were performed to (i) character- 
ize dose-efficacy responses of dipeptidyl peptidase 4 (DPP- 
4) inhibitors and TAK-875, and (ii) project and compare HbA1 c 
response with glimepiride, DPP-4 inhibitors (drug class), and 
TAK-875 up to 24 weeks. Model-based simulations provided a 
quantitative rationale for dose selection in phase 3. 



RESULTS 

Demographics and disposition 

This population PK analysis used 1,211 TAK-875 plasma 
samples collected from 286 patients with T2DM. Population 
PK-efficacy models were developed based on the data col- 
lected from all patients. Patients who received placebos were 
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Table 1 Summary statistics of baseline subject demographics and clinical laboratory measures 

Variable (continuous) Pharmacokinetics Pharmacokinetics-efficacy 





N 


Mean ± SD 


N 


Mean ± SD 


Age (years) 


286 


51.2± 10.5 


346 


51.5 ± 10.6 


BMI (kg/m 2 ) 


286 


31 .8 ±5.1 


346 


31.7±5.1 


Weight (kg) 


286 


86.3 ± 19.9 


346 


86.0 ± 20.0 


ALB, g/l 


286 


41 .7 ±3.2 


346 


41 .8 ±3.3 


ALP, U/l 


286 


84.8 ± 22.5 


346 


84.8 ± 24.4 


ALT, U/l 


286 


29.2 ± 16.7 


346 


29.4 ± 16.8 


AST, U/l 


286 


23.7 ± 11.4 


346 


24.0 ± 11.9 


CrCL, ml/min 


286 


130.9 ±43.3 


346 


129.7 ±44.0 


TBIL, umol/l 


286 


8.3 ±4.2 


346 


8.4 ±4.2 


DD, years 


285 


5.7 ±4.9 


344 


5.7 ±4.9 


GGT, U/l 


286 


40.7 ± 33.3 


346 


41 .4 ±33.4 


Metformin dose (mg) 


221 


1 ,696.2 ± 444.9 


273 


1,696.3 ±441.9 


BFPG (mg/dl) 


286 




346 


170.3 ±51.1 


Baseline HbA1c (%) 


286 




346 


8.4 ±0.93 


Variable (categorical) 


Category 


Percentage (%) 


Category 


Percentage (%) 


Gender 


Male 


47.2 


Male 


46.5 




Female 


52.8 


Female 


53.5 


Race 


African American 


9.4 


African American 


10.1 




Caucasian 


83.6 


Caucasian 


82.4 




Asian 


2.8 


Asian 


3.2 




Others 


4.2 


Others 


4.3 


Ethnicity 


Non-Hispanics 


30.8 


Non-Hispanics 


33.2 




Hispanics 


69.2 


Hispanics 


66.8 


Smokers 


Yes 


16.1 


Yes 


15.6 




No 


83.9 


No 


84.4 


Metformin 


Yes 


77.3 


Yes 


76.3 




No 


22.7 


No 


23.7 


Aspirin 


Yes 


11.9 


Yes 


13.0 




No 


88.1 


No 


87.0 


Angiotensin-converting enzyme Inhibitors 


Yes 


25.9 


Yes 


24.6 




No 


74.1 


No 


75.4 


Statins 


Yes 


15.0 


Yes 


16.5 




No 


85.0 


No 


83.5 



ALB, albumin; ALP, alkaline phosphate; ALT, alanine aminotransferase; AST, aspartate aminotransferase; CrCL, creatinine clearance; DD, disease duration; 
GGT, y-glutamyl transpepdidase; TBIL, total bilirubin. 



also included for PK-efficacy analysis. This population PK- 
efficacy analysis used 2,710 FPG and 1 ,381 HbA1 c samples 
collected from 346 patients with T2 DM. The overall population 
was predominantly Caucasian (82.4%), with mean age 52 
years (range 21-79) and mean body weight 86.0 kg (range 
49.5-172.7). Subject demographics and covariates are sum- 
marized in Table 1. 

PK analysis 

A two-compartment model (Supplementary Data online, 
control stream 3) with first-order absorption and elimination 
well described the exposure data (Figure 1). Interindividual 
variability (I IV) was estimated only for oral clearance (CL/F). 
Following covariate analysis, sex was identified as a statisti- 
cally significant covariate for TAK-875 CL/F as described in 
the following equation. 

CL/F = 0.752 + 0.308 x sex 

(sex= 1 for male and 0 for female) 



Males appear to have higher clearance as compared 
with females (Figure 2). Inclusion of sex as a covariate 
for CL/F of TAK-875 reduced the objective function value 
and I IV for CL/F in the final model by 14.52 and 7.6% 
coefficient of variation units, respectively. The parameter 
estimates along with their associated precisions (% rela- 
tive standard error for the final PK model for TAK-875) are 
presented in Table 2. 

Exposure-efficacy response analysis 

The effect of TAK-875 on FPG was well characterized by 
an indirect response model (Supplementary Data online, 
control stream 2) with stimulation of K out (first-order rate 
constant for elimination of glucose). I IV was estimated for 
the model predicted baseline (BL) and E max . The parameter 
estimates and their associated precision for the final FPG 
model for TAK-875 are presented in Table 2. Estimated E max 
of TAK-875 increased exponentially with observed baseline 
FPG (BFPG) whereas it increased linearly with increasing 
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aspartate aminotransferase (AST) levels as described in 
the following equation: 



E max = 0.366 x EXP (0.00746 x(BFPG- 163.5)) 
+ 0.00731x(AST-21) 



(2) 



The identified covariates (BFPG and AST) reduced the 
objective function value and I IV for E max in the final model by 
49.3 and 18.5% coefficient of variation units, respectively. 

The relationship between FPG and HbA1 c was implemented 
in the HbA1c model (Supplementary Data online, control 
stream 1 ) to describe the time course of HbA1 c (Supplemen- 
tary Figure S2 online). An exponential placebo model with 
factor for maximum placebo response (MPL) was included 



Dose 



Peripheral (V2) 
Compartment 



Gl 




Central (V1) 


depot 




Compartment 




Glucose 



HbA1c 



K A 



Figure 1 Population pharmacokinetics-efficacy model for TAK-875. 
Gl, gastrointestinal; HbA1c, glycosylated hemoglobin A1c. 

% Change in CL 
-100 -60 -20 20 60 100 



Gender 




in the model to empirically account for the effect of placebo 
on HbA1c levels observed in this short-term (12 weeks) trial. 
Due to the short trial duration, the half-life of MPL was fixed to 
the value (720 h) observed based upon the graphical analysis 
of data. IIV was estimated for the BL for HbA1c (BLA1), K MC , 
and MPL with a covariance term between MPL and K MC .The 
parameter estimates and their associated precision for the 
final HbA1c model for TAK-875 are presented in Table 2. 

The identified covariates (BFPG, sex, and disease dura- 
tion in years (DD)) reduced the IIV for BLA1, K MC , and MPL 
in the final model by 3.2, 29.8, and 7.9% coefficient of varia- 
tion units as compared with the base model without any 
covariates. BFPG was added as a covariate for the predicted 
HbA1c baseline during model refinement step. The relation- 
ship between identified covariates and efficacy parameters is 
presented in the equations below. 

BLA1 = (8.32 + 0.01 33x (DD - 4.61 )) x 
EXP (0.001 81 x (BFPG - 1 63.5)) 



(3) 



MPL = 0.0590 + 0.0363 x (sex) 



(4) 



The magnitude of the covariate effect on PK and efficacy 
parameters of TAK-875 are presented in Figure 2. Standard 
diagnostic plots for the developed population PK and effi- 
cacy models are shown in Supplementary Figure S3 online. 
Results of bootstrap evaluations and visual predictive check 
supported the robustness and predictive ability of the PK and 
efficacy models (Figure 3a,b, Table 2). 

Simulations outcome 

Simulation performed based on the developed population 
PK-efficacy model predicted mean reductions in HbA1c 

% Change in E max 
-120 -72 -24 24 72 120 



AST 



FPG 




262 



Covariate range represents 90% of the population 



% Change in maximum placebo response 
-65 -39 -13 13 39 65 



Sex 




Covariate range represents 90% of the population 



% Change in baseline HbA1c 
-65 -39 -13 13 39 65 



DD 



FPG 



0.39 I 14.75 



100 




262 



Covariate range represents 90% of the population 



Figure 2 Tornado Plot of covariate effects on pharmacokinetics and efficacy parameters of TAK-875. The reference line represents 
the typical subject in this population: a 52-year-old female subject with AST level of 21U/L, BFPG of 163.5 mg/dl and having diabetes 
for 4.61 years. Extreme values represent 5th and 95th percentile value observed for a covariate in this phase 2 study. AST, aspartate 
aminotransferase; BFPG, baseline fasting plasma glucose; DD, disease duration. 
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Table 2 Model and bootstrap parameter estimates with 95% bootstrap CI for 
the final TAK-875 PK, PK-FPG, and PK-efficacy model 



Parameter 


Model estimate 
(% RSE) 


Bootstrap 
estimate 


Bootstrap 

95% Cl a 


Pharmacokinetic analysis 








CL/F(l/h) 


0.75 (7.46) 


0.76 


0.67-0.98 


V1/F(l) 


5.86 (11.2) 


5.59 


3.858-72.2 


Q/F(l/h) 


0.833 


0.833 


NA 


V2/F(l) 


23.7 


23.7 


NA 


/((l/h) 


0.075 


0.075 


NA 


Gender on CL/F(l/h) 


0.31 (27.2) 


0.30 


0.14-0.55 


co 2 on CL/F 


0.488 (7.45) 


0.479 


0.276-1.22 


o 2 (exponential error) 


0.152 (3.67) 


0.149 


0.123-0.177 


Pharmacokinetic-FPG analysis 






BL (mg/dl) 


164(1.5) 


164 


159-168 




0.00542 (12.5) 


0.00543 


0.00419-0.0067 


^inax 


0.366 (16.2) 


0.366 


0.259-0.5006 


EC 50 (|jg/ml) 


3.16 (32.3) 


3.15 


1 .2-6.7 


BFPG on E max (mg/dl) 


0.00746 (13.1) 


0.00753 


0.00454-0.00971 


ASTonE max (U/l) 


0.00731 (44.2) 


0.00723 


0.00319-0.0122 


co 2 on BL 


0.0610 (9.3) 


0.0612 


0.0518-0.0721 


co 2 on E max 


0.0854 (53.9) 


0.0754 


0.00733-0.182 


o 2 (additive error) b 


0.0187(1.39) 


0.0186 


0.0163-0.0213 


Pharmacokinetic-efficacy analysis 






BLA1 (%) 


8.25 (0.5) 


8.26 


8.18-8.33 




0.00052 (1 1 .3) 


0.00053 


0.00042-0.00064 


HL(h) 


720 (fixed) 


NA 


NA 


MPL (NA) 


0.0590 (21 .7) 


0.0579 


0.0377-0.0794 


Gender on MPL 


0.0363 (32.1) 


0.0394 


0.0203-0.0631 


DD on BLA1 (%/year) 


0.0133 (55.6) 


0.0141 


0.0023-0.030 


BFPG on BLA1 (%/mg/dl) 


0.00181 (5.1) 


0.0018 


0.0016-0.0021 


co 2 on BLA1 


0.0057 (9.4) 


0.0057 


0.0047-0.0068 


co 2 on K MC 


0.95 (17.6) 


1.00 


0.658-1.55 


cov MPL ^ KA1c 


0.506 (35.0) 


0.486 


0.289-0.656 


co 2 on MPL 


0.305 (44.9) 


0.307 


0.168-0.540 


o 2 (additive error) b 


0.00164 (3.3) 


0.00156 


0.00131-0.00189 



Bootstrap convergence for PK, PK-PD, and PK-PD-efficacy were 99.9, 93.2, 
and 78%, respectively. 

AST, aspartate aminotransferase; BFPG, baseline FPG; BL, model predicted 
baseline; BLA1 , BL for HbA1c; CI, confidence interval; DD, disease duration; 
FPG, fasting plasma glucose; HbA1c, glycosylated hemoglobin A1c; HL, 
half-life of MPL; MPL, maximum placebo response; NA, not applicable; PD, 
pharmacodynamic; PK, pharmacokinetics; RSE, relative standard error. 
a CI calculated based on median values. b Additive error was used on 
log-transformed scale. 

levels from baseline following administration of 25 mg q.d. 
dosing as -0.94% (observed value = -0.84%) and -1.24% 
at the end of 3 and 6 months, respectively. Predicted reduc- 
tions from baseline HbA1c following 50 mg q.d. dosing were 
-1.16% (observed mean value = -1 .05%) and -1 .51% at the 
end of 3 and 6 months, respectively. 

Model-based meta-analysis 

With the model-based meta-analysis, time-response curves for 
DPP-4 inhibitors (therapeutic doses), glimepiride (therapeutic 
dose of 4 mg), and TAK-875 (selected therapeutic dose of 50 



mg) were projected up to 24 weeks (Figure 4a); and dose- 
efficacy response curves for DPP-4 inhibitors and TAK-875 
were characterized. Results indicated that therapeutic doses 
of DPP-4 inhibitors are associated with -80-90% of maximum 
effect on HbA1c, whereas 50 mg of TAK-875 is expected to 
produce -85% of maximum effect on HbA1c (Figure 4b). 

DISCUSSION 

One purpose of our analysis was to develop a population 
PK-efficacy model to describe the time course of TAK-875 
plasma concentrations and FPG and HbA1c levels follow- 
ing oral 6.25-200 mg q.d. doses in patients with T2DM for 
12 weeks. Furthermore, the main objective of the present 
analysis was to provide a quantitative rationale for selec- 
tion of doses for phase 3 trial using two pharmacometric 
approaches. 

The disposition kinetics of TAK-875 in patients with 
T2DM was best described using a PK model with first- 
order absorption and elimination processes. One- and 
two-compartment models were evaluated, with the latter 
PK model resulting in lower objective function value and 
Akaike Information Criterion values with substantially 
improved diagnostic plots. However, the two-compartment 
base PK model showed signs of instability during the 
bootstrap evaluation, which could indicate model over- 
parameterization. The model parameters describing the 
absorption and distribution of TAK-875 could not be pre- 
cisely estimated because of the sparse sampling in this 
phase 2 study. 10 It was also determined that the model was 
not sensitive to the particular values of intercompartmental 
clearance (Q/F), peripheral volume of distribution (V3/F), 
and first-order absorption rate constant (K a ). Therefore, 
these parameters were fixed to values estimated in an 
earlier population PK analysis performed using frequently 
collected PK samples in a multiple rising dose study in 
patients with T2DM. 9 Although potential for nonlineari- 
ties in the PK model (e.g., Michaelis-Menten elimination) 
was not tested, the lack of bias in the diagnostic plots 
suggested that a linear PK model adequately described 
the data and that TAK-875 exhibited linear PK over the 
evaluated dose range. Furthermore, approximate dose 
proportionality was previously demonstrated for TAK-875 
following oral administration of 25-800 mg single doses in 
healthy subjects 7 and 25-400 mg once-daily oral doses for 
1 4 days in patients with T2DM. 9 PK parameters of TAK-875 
were also constant over time after multiple administrations 
of 200 and 400 mg doses in healthy subjects. 8 These data 
together further support the appropriateness of consider- 
ing a linear PK model for TAK-875. 

Among baseline characteristics, only sex was identified 
as a statistically significant covariate on the CL/F of TAK- 
875, with -41% higher clearance in males (1.1 l/h) than in 
females (0.75 l/h). The estimated CL/F value for males in this 
study is consistent with the CL/F of 0.9-1 .7 l/h previously 
observed in healthy subjects. 7 TAK-875 primarily undergoes 
phase 2 metabolism via glucuronide conjugation. The results 
of this covariate analysis for PK of TAK-875 are not unex- 
pected, given that higher glucuronidation activity in males 
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Figure 3 (a) Visual predictive check for the fasting plasma glucose using final pharmacokinetics-efficacy model, (b) Visual predictive check 
for the HbA1c using final pharmacokinetics-efficacy model. FPG, fasting plasma glucose; HbA1c, glycosylated hemoglobin A1c. 



has previously been reported for drugs mainly metabolized 
by glucuronidation. 1112 On the basis of the simulations, males 
are expected to have maximum concentration and area under 
the curve values at steady state -25 and 29% lesser than 
those of females, respectively. This difference is expected to 
lead to <5 mg/dl and <0.01% differences in FPG and HbA1c 
levels between males and females over 12 weeks, and is not 
expected to have a clinically meaningful impact, given the 
dose response for efficacy and safety of TAK-875. 

Insulin-glucose homeostatic pharmacodynamic (PD) 
models reported in the literature 1314 were not evaluated 
as TAK-875 did not appear to have an apparent effect 
on fasting serum insulin based on the graphical exami- 
nation of fasting serum insulin levels vs. time obtained in 
this phase 2 study. This may be partially attributed to the 
glucose-dependent action of TAK-875. TAK-875 effects on 
glucose-stimulated insulin secretion leading to increases 
in the insulin levels are expected to become more notice- 
able in the state of elevated and rapidly changing glucose 
levels as seen postprandially.This will be explored in future 
studies using more frequent sampling to capture the early 
rise in insulin levels after intake of a standard meal and 
by developing a mechanistic model to account for com- 
plete insulin-glucose homeostasis. The placebo effect for 
FPG response was investigated as linear and exponential 



placebo response terms. However, the inclusion of placebo 
response did not improve model predictability and resulted 
in deviation from observed data in the diagnostic plots dur- 
ing PK-FPG analysis, because an inconsistent and highly 
variable placebo response was observed for FPG. There- 
fore, no placebo effect was included in the base and final 
PK-FPG models. 

The covariate analysis for TAK-875 PD parameters iden- 
tified observed BFPG and AST for E max as statistically sig- 
nificant covariates. The effect of BFPG on TAK-875 E is 

max 

consistent with the strong correlation between BFPG and 
E max observed for many antidiabetic drugs 15-17 ; the higher 
the baseline for FPG, the greater was the magnitude of 
TAK-875 response. On the basis of simulations over 12 
weeks, the changes in FPG and HbA1c levels from base- 
line were -8.5 mg/dl and -0.62%, respectively, for a subject 
with observed BFPG value of 100 mg/dl (5th percentile). 
By contrast, the changes in FPG and HbA1c levels from 
baseline were -62 mg/dl and -1.79%, respectively, with 
the observed BFPG value of 262 mg/dl (95th percentile). 
The effect of AST on E of TAK-875 was mainly driven 

max J 

by seven outliers with AST levels >70 U/l at baseline. On 
removal of these seven patients, baseline AST did not 
have statistically significant effect on E max of TAK-875 (P 
value > 0.005). 
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Figure 4 (a) Predicted mean changes in placebo adjusted HbA1 c up to 24 weeks for the 50 mg dose of TAK-875, the 4 mg dose of glimiper- 
ide and the therapeutic doses of DPP-4 inhibitors (linagliptin 5 mg, saxagliptin 2.5 mg, sitagliptin 100 mg, and vildagliptin 100 mg). Shaded 
area shows the range of mean HbA1c responses for the four DPP-4 inhibitors, (b) Dose-response curves for TAK-875 and DPP-4 inhibitors 
(reference drug class). Doses are shown as % of maximum clinically tested dose (i.e., highest clinically tested dose is 100%). Drug effect 
is shown as % of maximum drug effect on HbA1c (% E max ) at week 24. Shaded area shows range of % £ max associated with therapeutic 
doses. Vertical lines delimit range of therapeutic doses. Dashed lines delimit range of dose response curves for DPP-4 inhibitors. DPP-4, 
dipeptidyl peptidase 4; HbA1c, glycosylated hemoglobin A1c. 



Disease duration had no significant effect as a covariate 
on PK or efficacy parameters if controlled for BLA1 . However, 
this may be due in part to the limited range of disease dura- 
tion in the study population (0.39-14.75 years from diagno- 
sis). 18 No differences in efficacy parameters were observed 
for sex. This supports our interpretation that the difference 
found for PK between males and females is minor and not 
relevant to the PD effect of TAK-875. 

Of note, none of the comedications evaluated (metformin, 
statins, aspirin, and angiotensin-converting enzyme inhibi- 
tors) had a statistically or clinically significant effect on PK 
and efficacy parameters of TAK-875. This further indicates 
the limited potential for drug-drug interactions between TAK- 
875 and these drugs based on various TAK-875 interaction 
studies (data not shown). The lack of effect of metformin on 
the PD profile of TAK-875 can be attributed to the fact that 
patients were already on metformin treatment before enter- 
ing this study. 

On the basis of the simulation of HbA1c levels for 24 
weeks, the predicted placebo-adjusted mean reductions in 
HbA1c levels from baseline to the end of 6 months were 
-0.80 and -1.06% for 25- and 50-mg doses of TAK-875, 
respectively (Supplementary Figure S4 online). There- 
fore, TAK-875 50-mg q.d. dosing is expected to cause a 
decrease in HbA1c level from baseline as comparable with 
that of sulfonylureas (0.9-2%) 10 ' 19 - 23 and potentially greater 
than that of DPP-4 inhibitors. 24-28 On the basis of the simula- 
tions and observed data over 12 weeks, only minimal thera- 
peutic gain (reduction in HbA1c levels from baseline) can 
be expected at TAK-875 doses >50 mg. 10 Doses <25 mg 
may exhibit adequate pharmacological response; however, 
it is not expected to provide efficacy as comparable with 
the currently approved T2DM therapies (DPP-4 inhibitors, 
sulfonylureas, and glucagon like peptide-1 analogs). 19-22 



The following limitations should be considered for the 
empirical model presented, (i) The effect of TAK-875 on 
HbA1c levels did not reach plateau by month 3 in the pres- 
ent phase 2 study; therefore, the projection beyond 1 2 weeks 
may be slightly over or under estimated; (ii) a substantial 
proportion of subjects on placebo needed rescue medica- 
tion for glycemic control and were excluded from this analy- 
sis after the start of the rescue medications. As a result, 
the placebo response was not implemented in the PK-FPG 
model; (iii) the developed model assumes that BFPG and 
HbA1 c are at steady state at the time of randomization which 
may not be true slightly biasing overall simulated response; 
and (iv) a caution should be exercised for the projection on 
durability of TAK-875 response and extrapolation of placebo 
response beyond 12 weeks due to the short duration of cur- 
rent phase 2 trial. The model will be further optimized, if nec- 
essary, once long-term data are available from phase 3 trial 
to evaluate the accuracy of the prediction using the current 
model and assumptions considered during the development 
of the model. 

Applying a model-based meta-analysis approach, pro- 
jected HbA1c response with TAK-875 at the selected thera- 
peutic dose of 50 mg was consistent with that based on the 
exposure-efficacy response analysis. Results indicated that 
50 mg of TAK-875 is expected to produce a potentially greater 
HbA1c response than DPP-4 inhibitors and a response simi- 
lar to that of sulfonylureas after 24 weeks of treatment. Fur- 
thermore, characterized dose-response curves indicated 
that 50 mg of TAK-875 is expected to produce -85% of maxi- 
mum effect on HbA1c, which is consistent with values associ- 
ated with therapeutic doses of DPP-4 inhibitors (-80-90% of 
maximum effect on HbA1c). 

In conclusion, findings from pharmacometric analyses pro- 
vided a quantitative rationale for dose selection of 25 and 50 
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Figure 5 Schematics of study design. 

mg q.d. TAK-875. Both doses will be evaluated in the future 
clinical phase 3 program. 



METHODS 



Exposure-efficacy models. Population PK and efficacy analy- 
ses for TAK-875 were conducted using ADVAN4 TRANS4 
(PK) or ADVAN6TRANS1 (PK-efficacy) subroutines and first- 
order conditional estimation with r|-e interaction method as 
implemented in NONMEM VII. 29 Selection of base and final 
models was guided by the minimum objective function value, 
Akaike Information Criterion, visual inspection of diagnostic 
plots, and the precision of parameter estimates. 

One- and two-compartment models with first-order absorp- 
tion and linear elimination were evaluated as base structural 
PK models. Efficacy models were developed to characterize 
the time course of FPG and HbA1 c levels. A semimechanistic 
indirect response model. 3031 with stimulation of K out (first-or- 
der rate constant for removal of fasting glucose) as depicted 
in Figure 1 and described in differential equations below was 
used to describe changes in FPG following administration of 
TAK-875. 



Study design. Data used here were obtained from a phase 
2, randomized, double-blind, placebo- and active comparator- 
controlled, parallel-group, multicenter study. 10 The study was 
conducted to evaluate the efficacy, safety, and tolerability of 
TAK-875 q.d. treatment in patients with T2DM who were inad- 
equately controlled on a stable dose of metformin as mono- 
therapy for at least 8 weeks before screening, or who had been 
without chronic antidiabetic therapy within 8 weeks before 
screening with an 8-week documented diet and exercise plan. 
Eligible patients were randomized to seven treatment groups 
(Figure 5). The study was conducted in accordance with 
the guidelines on Good Clinical Practice and with the Ethi- 
cal Standard for Human Experimentation established by the 
Declaration of Helsinki. The protocol and amendments were 
reviewed and approved by the Institutional Review Board. 

Sample collection. Four blood samples (two trough and two 
nontrough, 4 ml each) were collected from each subject at two 
separate occasions at visit 7 (week 4) and visit 8 or 9 (week 6 
or 8) for quantitation of TAK-875 steady-state plasma concen- 
trations. TAK-875 plasma concentrations were measured using 
a validated ultraperformance liquid chromatography coupled 
with tandem mass spectrometry analysis. 7 FPG and HbA1c 
samples were collected from all patients at screening, base- 
line, day 7 (only FPG), day 14 (only FPG), and weeks 4, 6 (only 
FPG), 8, 10 (only FPG), and 12. FPG and HbA1c levels were 
measured using a standard laboratory method in a Clinical 
Laboratory Improvement Act (CLIA) certified central laboratory. 
Plasma pharmacokinetic, FPG, and HbA1c data collected from 
subjects on glimiperide were not included in the analysis. 

Pharmacokinetic software. NONMEM software version VII, 
level 1.0 (ICON Development Solutions, Ellicott City, MD) 
was used for nonlinear mixed-effect population PK-efficacy 
modeling. 29 NONMEM output was accessed through KIWI 
graphical interface (Cognigen Corporation, Buffalo, NY). 
Xpose version 4.0 (Uppsala University, Uppsala, Sweden) 
was used for generalized additive modeling analysis. Model- 
based meta-analyses and graphical plots were generated 
using S-PLUS version 8.0 (Tibco Software, Palo Alto, CA), 
Microsoft Excel 2007, and R version 2.12.1 (The R Founda- 
tion for Statistical Computing, Vienna, Austria). 



dFPG 
tit 



= Kin -Kout ><(1 + S)X FPG 



^ ^rnax 



xCp 



EC 50 +Cp 



BL: 



(5) 



where K m is the zero order rate constant for the produc- 
tion of FPG, E is the maximum effect of TAK-875 on the 

' max 

removal of FPG, EC 50 is TAK-875 concentration resulting in 
the 50% of maximum reduction in FPG levels from base- 
line, BL is model predicted baseline for FPG, and Cp is 
the concentration of TAK-875 following administration of 
different doses of TAK-875. Changes in HbA1c were mod- 
eled secondary to changes in FPG and the relationship 
between FPG and HbA1c was implemented in the PK-effi- 
cacy model 15 as depicted in Figure 1 and described in the 
equation below: 



dHbAlc 
df 



BLA1 = — ^ 

K Air 



= /C iG xFPG-HbA1cx/C A1c 



(6) 



K G is the first-order rate constant for the production of 
HbA1c, BLA1 is model predicted baseline for HbA1c, and 
K MC is the first-order rate constant for removal of HbA1c. 
On the basis of the graphical analysis, a consistent reduc- 
tion in HbA1c levels was observed in patients who were on 
placebo. As a result, a placebo factor (PLAC), as defined by 
the equation below, was included in the differential equation 
to account for the effect of lifestyle intervention or placebo on 
HbA1c levels. 



PLAC=1-MPLx(l-e' 



((-Ln(2)/HL)xtime) 



(7) 



dHbAlc 
tit 



K iG x (PLAC) x FPG - HbA1 c x K MC (8) 



HL is the half-life of MPL, the maximum placebo response, 
and PLAC refers to overall placebo effect. I IV associated 
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with PK and efficacy parameters of TAK-875 was modeled 
assuming a log-normal distribution. 



g(x) 



.(f)- Dose. 



Dose +ED cn ^ 

// 50, drug 



(11) 



:P pop XeX P(^) 



(9) 



where P is the estimated parameter value for individual 
i, P pop represents the typical population estimate for the 
parameter, and r| i is the deviation of P from P pop . A diagonal 
covariance matrix was implemented for PK and PK-FPG anal- 
yses whereas a block matrix was implemented for the HbA1c 
analysis. 

Residual variability was modeled using a proportional error 
term for PK (non log-transformed TAK-875 plasma concentra- 
tion) and a log-normal error term (additive on log-transformed 
scale) for PK-FPG and HbA1c analyses. Only clinically rel- 
evant and biologically plausible covariate-parameter rela- 
tionships were explored in the covariate analysis. In addition, 
comedications taken by more than 1 0% of the patient popula- 
tion during 80% of the trial period were also evaluated in the 
covariate analysis. Potential covariates identified by general- 
ized additive modeling and graphical analysis (plot of r\. vs. 
covariate) were then tested using stepwise forward addition 
(P < 0.05) followed by stepwise backward elimination proce- 
dure (P< 0.005). 32 

Model qualification was conducted using standard non- 
parametric bootstrap (n = 1 ,000) and visual predictive check 
to evaluate the precision of model parameter estimates and 
robustness of the final PK and efficacy models. 33 



Model-based meta-analysis. A model-based meta-analysis 
used controlled clinical trials from the medical literature, the 
Food and Drug Administration, and the European Medicines 
Agency. 34 Clinical outcomes data from 74 prospective ran- 
domized clinical trials, with >36,300 patients and 208 different 
randomized treatment arms, were included in this model- 
based meta-analysis on efficacy of glimepiride and DPP-4 
inhibitors (alogliptin, linagliptin, saxagliptin, sitagliptin, and 
vildagliptin). For TAK-875, data from a randomized placebo- 
controlled phase 2 study (excluding glimepiride arm) were 
evaluated. 10 

Relationships between baseline HbA1c and doses of 
DPP-4 inhibitors and TAK-875 were analyzed using a non- 
linear regression implemented in the Generalized Nonlinear 
Least Squares routines of R version 2. 12.1. 34-36 Changes in 
HbA1c in a given treatment arm (/) of a trial (/) (HbA1c.) at 
time T were modeled as a function of the placebo response 
for HbA1c in that trial (E QjT ), i.e., the intercept, and a dose- 
response relationship for treatment effects on HbA1c (g(x)). 
Covariates (X) and trial-specific model parameters (0.) were 
included and evaluated according to the following general 
structure: 

HbA1 c ijT = E 0JJ + g(Drug v Dose v X v 0 n T) (10) 

where E 0/T is a nonparametric placebo (i.e., reference) effect, 
which takes a different value for each study-time combination 
for each trial. 35-37 Dose-response relationships for random- 
ized treatments were characterized as follows: 



where E maxciass is the maximal drug effect, reflecting the maxi- 
mal reduction in placebo-adjusted HbA1c (i.e., difference in 
response between placebo and active treatment) at a given 
time (t). 35-37 A different E max was estimated for DPP-4 inhibi- 
tors (class effect), glimepiride, and TAK-875. This assump- 
tion was tested by allowing for a different E max for each DPP-4 
inhibitor. Dose is the total daily dose and ED 50 is the dose to 
achieve 50% of E .A different ED Kn was estimated for each 

max 50 

drug. 

Factors such as background treatment and base- 
line HbA1c were included as covariates on the param- 
eter E in the model. Additional random between-trial 

max 

heterogeneity in the relative effect between arms was 
accounted for by the trial-specific model parameters that 
were assumed to be normally or log-normally distributed 
with between-trial variance Q. Model selection was based 
on a Log likelihood ratio test and the confidence inter- 
val of the parameter estimate at an acceptance P value 
of 0.01. Confidence intervals of the parameter estimates 
were derived from the variance matrix of the parameter 
estimates. 3536 

For DPP-4 inhibitors, both apparent onset and durability 
of drug-related effects on HbA1c were estimated. For TAK- 
875 (due to lack of long-term data), only the onset effect on 
HbA1c was estimated. To take a conservative approach, the 
"offset" of TAK-875 related effects (i.e., the durability of TAK- 
875-related effects) on HbA1c was assumed to be the same 
as that of DPP-4 inhibitors. 37 

Simulations for projection and comparison of long-term effi- 
cacy. The final PK and efficacy models were used to perform 
simulations to predict TAK-875, FPG, and HbA1c concentra- 
tions following administration of 25- and 50-mg q.d. doses of 
TAK-875 over a period of 6 months (24 weeks). The simula- 
tions were used to support dose selection for future studies 
and also to compare the efficacy of TAK-875 with currently 
approved pharmacological T2DM therapies. The simulations 
were performed using a population of 1 ,000 patients receiv- 
ing 25- and 50-mg multiple q.d. doses of TAK-875. These 
sets of simulations assumed durability of efficacy responses 
beyond 12 weeks, which was the duration of this clinical trial 
used to generate data utilized to develop PK and efficacy 
models for TAK-875. 

The developed meta-analysis model for HbA1c utilizing 
publicly available clinical trial outcomes data was also used 
to project time-response curves for DPP-4 inhibitors (refer- 
ence drug class) and TAK-875 at the selected therapeutic 
dose of 50 mg up to 24 weeks. Baseline HbA1c was set to 
8.5%, i.e., the baseline value observed in the phase 2 study. 10 
A time-response curve for DPP-4 inhibitors (drug class) at 
therapeutic doses was projected up to 24 weeks (linagliptin 
5 mg, saxagliptin 2.5 mg, sitagliptin 100 mg, and vildaglip- 
tin 100 mg). In addition, dose-efficacy response curves of 
DPP-4 inhibitors (reference drug class) and TAK-875 were 
also characterized. 



CPT: Pharmacometrics & Systems Pharmacology 



GPR40 agonist, TAK-875, Population PK and Efficacy 

Naik eta\ 



Acknowledgments. Funding for this study was provided by 
Takeda Pharmaceuticals North America, Inc. This study was 
supported by Takeda Global Research & Development Center, Inc. 
Deerfield, IL. 

Author contributions. M.V., H.N., C.C., M.R, and E.L. wrote the 
manuscript; M.V., H.N., and E.L. designed the research; M.V., H.N., 
J.L., and E.L. performed the research; M.V., H.N., J.L., C.C., M.R, 
and E.L. analyzed the data. 

Conflict of interest. H.N., E.L., C.C., and M.V. are employees of 
Takeda Global Research & Development Center, Inc., Deerfield, IL; 
J.L. and M.R are employees of Quantitative Solutions, Inc., Bridge- 
water, NJ. 



1. Wild, S., Roglic, G., Green, A., Sicree, R. & King, H. Global prevalence of diabetes: esti- 
mates for the year 2000 and projections for 2030. Diabetes Care 27, 1047-1053 (2004). 

2. National Diabetes Statistics, 2007. National Institute of Health, National Institute of Diabe- 
tes and Digestive and Kidney Diseases (NIDDK). Publication No. 08-3892. 2008. 

3. American Diabetic Association. National Diabetes Fact Sheet, 201 1 . <http://www.diabetes. 
org/in-my-community/local-offices/miami-florida/assets/files/national-diabetes-fact-sheet. 
pdf>.(201 1). Accessed 1 March 2012. 

4. Saydah, S.H., Fradkin, J. & Cowie, C.C. Poor control of risk factors for vascular disease 
among adults with previously diagnosed diabetes. JAMA 291, 335-342 (2004). 

5. Negoro N. et al. Discovery of TAK-875: a potent, selective, and orally bioavailable GPR40 
agonist. ACS Med. Chem. Lett. 1 , 290-294 (2010). 

6. Tsujihata, Y. et al. TAK-875, an orally available G protein-coupled receptor 40/free fatty 
acid receptor 1 agonist, enhances glucose-dependent insulin secretion and improves both 
postprandial and fasting hyperglycemia in type 2 diabetic rats. J. Pharmacol. Exp. Ther. 
339, 228-237 (2011). 



7. Naik, H. et al. Safety, tolerability, pharmacokinetics, and pharmacodynamic 
properties of the GPR40 agonist TAK-875: results from a double-blind, placebo- 
controlled single oral dose rising study in healthy volunteers. J. Clin. Pharmacol. 52, 
1007-1016 (2012). 

8. Matsuno K., Hirayama M., Araki T., Dote N., Kondo T. & Nakamura K. Pharmacokinetics, 
safety and tolerability of single and multiple doses of TAK-875, a novel GPR40 agonist, in 
Japanese healthy male subjects. Diabetes 59, A193 (2010). 

9. Leifke, E. et al. A multiple-ascending-dose study to evaluate safety, pharmacokinetics, 
and pharmacodynamics of a novel GPR40 agonist, TAK-875, in subjects with type 2 dia- 
betes. Clin. Pharmacol. Ther. 92, 29-39 (2012). 

10. Burant, C.F. et al. TAK-875 versus placebo or glimepiride in type 2 diabetes mellitus: 
a phase 2, randomised, double-blind, placebo-controlled trial. Lancet 379, 1403-1411 
(2012). 

1 1 . Tanaka, E. Gender-related differences in pharmacokinetics and their clinical significance. 
J. Clin. Pharm. Ther. 24, 339-346 (1999). 

12. Greenblatt, D.J. & von Moltke, L.L. Gender has a small but statistically significant 
effect on clearance of CYP3A substrate drugs. J. Clin. Pharmacol. 48, 1350-1355 
(2008). 

13. Ribbing, J., Hamren, B., Svensson, M.K. & Karlsson, M.O. A model for glucose, insulin, 
and beta-cell dynamics in subjects with insulin resistance and patients with type 2 diabetes. 
J. Clin. Pharmacol. 50, 861-872 (2010). 

14. de Winter, W. et al. A mechanism-based disease progression model for comparison 
of long-term effects of pioglitazone, metformin and gliclazide on disease processes 
underlying Type 2 Diabetes Mellitus. J. Pharmacokinet. Pharmacodyn. 33, 313-343 
(2006). 

1 5. Rohatagi, S. et al. Model-based development of a PPARgamma agonist, rivoglitazone, to 
aid dose selection and optimize clinical trial designs. J. Clin. Pharmacol. 48, 1420-1429 
(2008). 

1 6. Marino, M.T. et al. Pharmacokinetics and pharmacodynamics of inhaled GLP-1 (MKC253): 
proof-of-concept studies in healthy normal volunteers and in patients with type 2 diabetes. 
Clin. Pharmacol. Ther. 88, 243-250 (2010). 

17. Frey, N., Laveille, C, Paraire, M., Francillard, M., Holford, N.H. & Jochemsen, 
R. Population PKPD modelling of the long-term hypoglycaemic effect of gliclazide given 
as a once-a-day modified release (MR) formulation. Br. J. Clin. Pharmacol. 55, 147-157 
(2003). 

18. Verma, M., Paneri, S., Badi, P. & Raman, P.G. Effect of increasing duration of diabetes 
mellitus type 2 on glycated hemoglobin and insulin sensitivity. Indian J. Clin. Biochem. 21, 
142-146 (2006). 

1 9. Rodbard, H.W. et al. American Association of Clinical Endocrinologists medical guidelines 
for clinical practice for the management of diabetes mellitus. Endocr. Pract. 13(suppl. 1), 
1-68 (2007). 

20. Nathan, D.M. et al. Medical management of hyperglycemia in type 2 diabetes: a 
consensus algorithm for the initiation and adjustment of therapy: a consensus statement 
of the American Diabetes Association and the European Association for the Study of 
Diabetes. Diabetes Care 32, 193-203 (2009). 

21. Schade, D.S., Jovanovic, L. & Schneider, J. A placebo-controlled, randomized study of 
glimepiride in patients with type 2 diabetes mellitus for whom diet therapy is unsuccessful. 
J. Clin. Pharmacol. 38, 636-641 (1998). 

22. Simonson, D.C., Kourides, I.A., Feinglos, M., Shamoon, H. & Fischette, C.T. Efficacy, 
safety, and dose-response characteristics of glipizide gastrointestinal therapeutic system 
on glycemic control and insulin secretion in NIDDM. Results of two multicenter, random- 
ized, placebo-controlled clinical trials. The Glipizide Gastrointestinal Therapeutic System 
Study Group. Diabetes Care 20, 597-606 (1997). 

23. Amaryl [package insert]. Bridgewater, NJ; Sanofi-Aventis, 2012. 

24. Raz, I., Hanefeld, M., Xu, L, Caria, C, Williams-Herman, D. & Khatami, H. Efficacy and 
safety of the dipeptidyl peptidase-4 inhibitor sitagliptin as monotherapy in patients with type 
2 diabetes mellitus. Diabetologia 49, 2564-2571 (2006). 

25. Rosenstock, J., Aguilar-Salinas, C, Klein, E., Nepal, S., List, J. & Chen, R. Effect of saxa- 
gliptin monotherapy in treatment-naive patients with type 2 diabetes. Curr. Med. Res. Opin. 
25, 2401-2411 (2009). 

26. Lukashevich, V., Schweizer, A., Shao, Q., Groop, P.H. & Kothny, W. Safety and efficacy of 
vildagliptin versus placebo in patients with type 2 diabetes and moderate or severe renal 
impairment: a prospective 24-week randomized placebo-controlled trial. Diabetes. Obes. 
MetabAZ, 947-954 (2011). 

27. Onglyza [package insert], Prinston, NJ; Bristol-Myers Squibb Company, 201 1 . 

28. Januvia [package insert], Whitehouse Station, NJ; Merck & Co, 2012. 

29. Nonmem. NONMEM Users Guide: Introduction to NONMEM 7. (ICON Develop- 
ment Solutions, Ellicott City, MD, 2010). 

30. Landersdorfer, C.B. & Jusko, W.J. Pharmacokinetic/pharmacodynamic modelling in diabe- 
tes mellitus. Clin. Pharmacokinet. 47, 417-448 (2008). 

31 . Lin, S. & Chien, Y.W. Pharmacokinetic-pharmacodynamic modelling of insulin: comparison 
of indirect pharmacodynamic response with effect-compartment link models. J. Pharm. 
Pharmacol. 54, 791-800 (2002). 

32. Wahlby, U., Jonsson, E.N. & Karlsson, M.O. Assessment of actual significance levels 
for covariate effects in NONMEM. J. Pharmacokinet. Pharmacodyn. 28, 231-252 
(2001). 



Study Highlights 

WHAT ISTHE CURRENT KNOWLEDGE ONTHE 
TOPIC? 

GPR40 or free fatty acid receptor 1 , expressed in pancre- 
atic p, has become a drug target for treatment of T2DM. 
GPR40 agonists increase glucose-stimulated insulin 
secretion and have demonstrated beneficial effects in 
various preclinical models of T2DM. TAK-875 is a novel 
and first GPR40 agonist in phase 3 clinical development. 
TAK-875 produces a similar glucose-lowering effect as 
sulphonylureas with a considerably lower potential for 
hypoglycemia. 

WHAT QUESTION DIDTHIS STUDY ADDRESS? 

The analysis describes the pharmacometric approaches 
that guided the selection of TAK-875 dose for phase 3 
studies in patients with type 2 diabetes. 

WHATTHIS STUDY ADDSTO OUR KNOWLEDGE 

This analysis demonstrates the importance of pharmaco- 
metric approaches for rationale therapeutic dose selection 
for phase 3 trial. 

HOWTHIS MIGHT CHANGE CLINICAL 
PHARMACOLOGY ANDTHERAPEUTICS? 

The analysis presented emphasizes the importance of 
model-based drug development approach for the clinical 
development of drugs. The pharmacometric approaches 
implemented in this analysis provided a quantitative ratio- 
nale for the dose selection and a comparative analysis to 
compare the efficacy of TAK-875 with marketed product 
(DPP4i). 
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